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A NUMERICAL INTEGRATION SCHEME TO DETERMINE 
HEMISPHERIC EMITTANCE, SOLAR ABSORPTANCE 
AND EARTH INFRARED ABSORPTANCE FROM 
SPECTRAL REFLECTANCE DATA 


INTRODUCTION 


The total hemispheric emittance (€ ) , solar absorotance {a ) , and 

X s 

earth infrared absorptance (a ) of a selected coating are three of the 

IH 

quantities required to completely specify the optical condition of a surface that 
is to be used in a space application. This paper describes a numerical method 
of computing these required quantities from experimentally measured values 
for the spectral reflectance of a coating. 

The result of this method is a numerical integration scheme to calculate 
€ , a , and a: that can easiJy be programmed on a computer. The inputs 

X S XX V 

to the computer program are the spectral reflectance and transmittance or the 
spectral absorptance and transmittance. The output of this program for 
several thermal control coatings is given in Appendix A. 


TOTAL EMITTANCE 


The total emittance of a surface is defined as 


e 


T 


Energy total emitted by surface at temperature, T 
Energy total emitted by a black body at same terriperature 


thus, 



1 1 ) 


where 


a 


A 


spectral absorptance of the surface, 


A) = spectral energy distribution emitted by a black body at 
temperature T, and 


A = wavelength. 

Since the input data to this equation are experimental, the integration 
must be numerical, so the equation for total emittance becomes 


€ 


T 



y o S (A) AA . 
A=0 ^ ^ 


( 2 ) 


where E™, = f S (A)dA = total black body energy at temperature T. 

The method of ’’weighted ordinates" to calculate the total emittance 
was described by Atkins (1). In this method, equal values of A A are used in 
equation ( 2) . 

The method described here is the "100 selected ordinates" which 
reduces equation (2) to a simpler form. 

If the O' are picked such that each increment A A contains one percent 

A 

of the total black body energy, E^, then the sum becomes 

00 OO JT 

1 - ' 3 . 

A — U A — 0 


9 




but, since is independent of A 


V a E (A)AA 
A = 0 ^ 




Equa. 4 


The total emittance equation now becomes 


1 V 

= ' a 

^ A=0 ^ 


There is an a for each E • therefore, there are 100 indep^mdent ot ’s. 
Thus, A A A 


100 
y a 

100 m 

m=l 


The equation for total emittance has now been reduced to a sun\ of the 
spectral absorptances taken in a special way. 

Since the absorptivity is spectrally dependent, the wavelength of the 

a 's must be found, 
m 

In the typical black body curve seen in Figure 1, one percent of the 
energy is between the lines A and C. However, the that will be selected 

to represent this energy segment will not be at A. or A , but, at the mid- 

IK 

point, A or at the 0.5 percent interval. Therefore, the desired wavelengths 
will fall at 0. 5 percent, 1. 5 percent, etc. 


Black body curves vary with the temperature, but they are similar in 


that 


A. T = k. 
1 1 


(i = 1, 2, 3 ....) , 




ABC 




FIGURE 1. TYPICAL BLACK BODY CURVE 


where 

= wavelength 

T = black body temperature 

k. = constant 

1 


The k.'s for the 100 energy points 0. 5 percent, 1. 5 percent, 2. 5 percent, 

etc. , have been computed and are listed in Appendix B. Thus, for a given black 
body temperature, the wavelengths for the desired absorptances may be easily 
calculated. 

The experimental points for absorptance will seldom fall on the wave- 
lengths given by equation (7) ; therefore, a linear interpolation between the data 

points is used to find the a 's. 

m 


The spectral absorplance data that can be obtained with the equipment 
used does not cover the entire range of the black btxly curve, so the linear 
interpolation scheme did not give all the 

Three techniques were used to obtain the beyond the endpoints. 

The first of these is called ”no extrapolation. " In this method the limits on 
the integratioi; in equation (1) are reduced to the range covered by the experi- 
mental data. 


Then, for this case, equation ( 11 


is changed to 


A2 

y S_tAl AA 
At 

_ A=Ai ^ 1^ 

^T A^ lOOA 

^ S^lAl A A 
A=Ai 



f 


where Aj and A 2 are the endpoints of the data and 
A 2 



In the output of this program A }'epresents the percent of the black body 
curve covered by the data. The complement of A (in percent) is printed by 
tiie computer program along with the enr.ittance values under the heading 
"percent not covered. " 

In the second method, called "straight extrapolation," the value of a 
beyond the endpoints is taken to be the vtJue of at the endpoint ( Fig. 2; . 

The last method, called the "curved extrapolation," fits a parabola to 

the last three points at each end of the data ( Fig. 3) and uses this curve for 

choosing a in this region. In this method, the value of a is forced to 
m m 

stay between 0 and 1 since these values are the minimum and maximum for 

absorptance. 


ax 





f 


FIGURE 2, ’’STRAIGHT EXTRAPOLATION” SCHEME 

When all the are found and sumnied from K = 1 to 100 in equation 

(6) for all three endpoint cases, the results obtained will be three values of 
the total errittance. A comparison of the three values will yield a range of 
values for but the ”no extrapolation” case will usually provide a more 

useful value for the total emittance. 


SOLAR ABSORPTANCE 


The total solar absorptance of a surface is calculated the same way 
as the total emittance, since 


(r 


s 



s 


f o S A) dX « — 
J \ a E 

0 s 


oo 


I 

>-0 


Of. S (A) AA . 
A s 


6 


■i 


Qx 



FIGURE 3. "CURVED EXTRAPOLATION" SCHEME 

where 

or^ = total solar absorptance 

rj; = spectral absorptance of the sample 
A 

S ( A) = spectral solar energy distribution 
s 

E = integrated solar enerj^y 
s 

'ITius, if Jonnson's solar irradiance curve (2J ( Fig. 4) is divided into 100 
equal energy bands (Appendix B) , as in the total emittance case, and the 
Qf^^'s determined as before, the total solar absorptance is 

^ 100 

(X = — ' V (X 

s 100 m 

m-1 
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FIGURE 4. JOHNSON'S SOLAR UGiADIANCE CURVE [2] 

EARTH INFRARED ABSORPTANCE 


The earth absorptance, as in Uie solar ease, is defined by the 
equation 

OO OC 

o; = — “ f a ^ dA ^ ~~ V c* S 

m E J \ IR^^ E A IR 

0 

where 

Q! = spectral absorptance 

A, 

S^ = spectral earth IR emittance distribution 


E 


integrated earth IR emittance 



I 

7 If the energy under Johnson’s earth infrared emission curve ( Fig. 5 ) 

^ is divided into one percent energy segments (Appendix B) , and an absorpiance 

> Of ^ determined for each one percent energy segment, then tlie expression for 

f total earth infrared absorptance becomes 

f-' 

A ' .100 



•«>AVE LENGTH (^) 

FIGURE 5. JOHNSON'S EARTH INFIIARED EMISSION CURVE [2] 

APPLICATION TO SPECTRAL DATA 


The spectral data used as input to these equations are taken on two 
different spectrometers. The first is an integrating sphere which provides 
total hemispherical reflectance measurements from 0.295 microns to 2.0 
microns. This range covers all but four percent of the solar spectrum; 






therefore, integrating sphere measurements are sufficient to calculate solar 
absorptance, but not total emittance or earth n< absorptance. These require 
measurements in the longer wavelength region; therefore, a second spectrom- 
eter, a heated Ho hi raum (cavity), is used to make total hemispherical measure 
ments in the range from 2. 0 microns to 25 microns. This range is sufficient 
to cover Uie majority of the energy under the earth infrared curve and the 
black body curves (black body temperatures from approximately 200 to 500® K) . 

CONCLUSION 

Therefore, for all three cases (€ , a , a ) , the equations reduce 

. X S XX V 

to 

^ 100 

^'^i 100 ^ ^m ’ 

m=l 

where X. = o or with the a and AA. chosen in a special way. 

1 T s lit m t' j 

This enables easy programming of the problem on a computer with the required 

constants to find the A^^’s stored as permanent data. This program will then 

require only spectral data and black body temperatures to compute e , o; , 

X s 
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APPENDIX B 


ONE-PERuENT ENERGY BANDS 


TABLE 1. ONE PERCENT ENERGY BANDS FOR 
JOHNSON'S SOLAR CURVE 


%E 

\ 

%E 


. 5 

. 2793 

30. 5 

. 5500 

1.5 

. 3053 

31. 5 

.5575 

2.5 

.3235 

32.5 

. 5650 

3.5 

. 3355 

33. 5 

. 5721 

4.5 

. 3484 

34.5 

.5800 

5.5 

.3603 

35.5 

. 5871 

6.5 

.3714 

36. 5 

. 5950 

7.5 

.3824 

37.5 

. 6025 

8.5 

.3945 

38.5 

.6108 

9. 5 

.4038 

39. 5 

.6185 

10. 5 

.4110 

40. 5 

. 6267 

11. 5 

.4186 

41. 5 

. 6305 

12. 5 

.4258 

42. 5 

.6433 

13. 5 

.4336 

43.5 

.6417 

14. 5 

.4407 

44.5 

. 6600 

15. 5 

.4475 

45.5 

. 6691 

16. 5 

.4573 

46. 5 

. 6782 

17. 5 

.4600 

47. 5 

. 6873 

18. 5 

. 4669 

48.5 

. 6970 

19. 5 

.4731 

49.5 

. 7070 

20. 5 

.4794 

50.5 

.7170 

21.5 

. 4864 

51.5 

.7270 

22. 5 

.4931 

52.5 

.7378 

23.5 

. 5000 

0^ . 5 

.7980 

24. 5 

. 5071 

54.5 

.7595 

25.5 

. 5143 

55.5 

.7714 

26. 5 

. 5217 

56.5 

. 7833 

27. 5 

. 5293 

57. 5 

. 7952 

28. 5 

. 5364 

58. 5 

. 8079 

29. 5 

. 5431 

59. 5 

. 8210 


TABLE 1 (Concluded) 


%E 


%E 

2 _ 

GO. 5 

. 8342 

80. 5 

1. 1974 

Gl. 5 

. 8474 

81. 5 

1. 2273 

G2. 5 

. 8618 

82. 5 

1.2576 

G3. 5 

. 8765 

83. 5 

1. 2879 

G4. 5 

. 8912 

84. 5 

1.3231 

G5. 5 

. 9100 

85. 5 

1.3615 

66. 5 

. 9233 

86. 5 

1.4000 

G7. 5 

. 9400 

87.5 

1.5476 

G8. 5 

. 9571 

88. 5 

1.5952 

60. 5 

. 9750 

89. 5 

1. 6500 

70. 5 

. 9929 

90.5 

1.7070 

71. 5 

1.0125 

91.5 

1.7769 

72. 5 

1. 0333 

92.5 

1. 8535 

73. 5 

1.0542 

93. 5 

1. 9565 

74. 5 

1.0750 

94.5 

2. 0887 

75. 5 

1.0958 

95.5 

2.2588 

76. 5 

1. 1205 

96. 5 

2.4820 

77. 5 

1. 1461 

97.5 

2. 8136 

78. 5 

1. 1718 

98. 5 

2.3846 

79. 5 

1. 1718 

99. 5 

4. 9667 


%E 

TABLE IL ONE PERCENT ENERGY BANDS FOR 
JOHNSON’S EARTH IR CURVE 

%E 

. 5 

4. 92700 

11.5 

8. 83400 

1.5 

6. 53800 

12. 5 

8. 92600 

2. 5 

7. 61900 

13. 5 

9. 02000 

3. 5 

8.02300 

14.5 

9. 11400 

4. 5 

8. 16400 

15.5 

9.21000 

5.5 

8.27100 

16. 5 

9.31700 

G. 5 

8. 37000 

17.5 

9. 44300 

7. 5 

8.46600 

18. 5 

9. 60300 

8. 5 

8. 55800 

19. 5 

9. 85600 

9. 5 

8.65100 

20.5 

10.01100 

10. 5 

8.74300 

21. 5 

10. 11200 
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TABLE II (Continued) 


%E 

iL 

%E 

X 

i_ 

22.5 

10.20200 

54 . 5 

14.29800 

23.5 

10.29200 

55.5 

14.57200 

24.5 

10.38100 

56.5 

14.92500 

25.5 

10.4G900 

57.5 

15.28200 

2G.5 

10.55800 

58.5 

15.64500 

27.5 

10.64G00 

59.5 

16.00900 

28.5 

10.73500 

60 . 6 

16.34400 

29 . 5 

10.82300 

61.5 

16.65500 

30.5 

10.91100 

62.5 

16.95700 

31.5 

11.00000 

63.5 

17.24800 

32.5 

11.08800 

64.5 

17.51700 

33.5 

11. 17700 

65.5 

17.76800 

34.5 

11.2G700 

66.5 

18.01500 

35.5 

11.35700 

67.5 

18.27800 

3G.5 

11.44700 

68.5 

18. 58200 

37.5 

11.53800 

69 . 5 

18.90500 

38.5 

11.62900 

70.5 

19. 24700 

39.5 

11.72400 

71.5 

19. 60800 

40.5 

11.82200 

72.5 

19.98900 

41.5 

11.92800 

73.5 

20.39600 

42.5 

12.04700 

74.5 

20.83300 

43.5 

12.20700 

75.5 

21.30900 

44.5 

12.41100 

76.5 

21. 83600 

45.5 

12.61800 

77.5 

22.41500 

4G.5 

12.82100 

78.5 

23. 05000 

47.5 

13.01100 

79.5 

23. 72400 

48.5 

13. 19100 

80.5 

24.44000 

49.5 

13.36600 

81.5 

25. 20500 

50. 5 

13.54000 

82.5 

26. 02900 

51.5 

13.71600 

83.5 

27. 81400 

52.5 

13.89800 

84.5 

28.78300 

53.5 

14. 09000 




For a black body emission curve of temperature T 


X. T = K. . 
1 1 
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The following table is a listing of the K^’s at the 0.5%, 1. 5%, 2.5% 
99. 5% points on the curve making it possible to determine the X.'s at the 
desired energy points for any teniperature T. 


TAbLE UL BLACK BODY CONSTANTS 


% 


% 


. 5 

1322. 

31.5 

3186. 

1.5 

1534. 

32.5 

3231. 

2.5 

1662. 

33.5 

3277. 

3.5 

1762. 

34.5 

3323. 

4.5 

1846. 

35.5 

3369. 

5.5 

1920. 

36.5 

3415. 

6.5 

1989. 

37.5 

3462. 

7.5 

2050. 

38.5 

3510. 

8.5 

2111. 

39.5 

3558. 

9.5 

2168. 

40.5 

3607. 

10.5 

2222. 

41.5 

3656. 

11.5 

2274. 

42.5 

3706. 

12.5 

2325. 

43.5 

3757. 

13.5 

2374. 

44.5 

3809. 

14.5 

2423. 

45.5 

3861. 

15.5 

2470. 

46.0 

3914. 

16.5 

2517. 

47.5 

3968. 

17.5 

2563. 

48.5 

4023. 

18.5 

2609. 

49.5 

4079. 

19.5 

2654. 

50. 5 

4136. 

20.5 

2698. 

51.5 

4194. 

21.5 

2743. 

52.5 

4254. 

22.5 

2787. 

53.5 

4314. 

23.5 

2831. 

54.5 

4377. 

24.5 

2876. 

55.5 

4440. 

25.5 

2920. 

56.5 

4505. 

26.5 

2964. 

57.5 

4572. 

27.5 

3008. 

58.5 

4640. 

28.5 

3052. 

59.5 

4710. 

29.5 

3097. 

60.5 

4782. 

30.5 

3141. 

61.5 

4856. 
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TABLE III 

% 

G2.5 

63.5 

64.5 

65.5 

66.5 

67.5 

68.5 

69.5 

70.5 

71.5 

72.5 

73.5 

74.5 

75.5 

76.5 

77.5 

78.5 

79.5 

80.5 

81.5 

82.5 

83.5 

84.5 

85.5 

86.5 

87.5 

88.5 

89.5 

90.5 

91.5 

92.5 

93.5 

94.5 

95.5 

96.5 

97.5 

98.5 

99.5 


( Concluded) 

ilL 

4932. 
5010. 
509 1 . 
5175. 
5262. 
5351. 
5444. 
5541. 
5641. 
5745. 
5854. 
5968. 
6087. 
6212. 
6343. 
6482. 
6628. 
6783. 
6948. 
7123. 
7311. 
7514. 
7732, 
7969. 
8228. 
8513. 
8829. 
9183. 
9153. 
10042 . 
10577. 
11215. 
11996. 
12990. 
14327. 
19724. 
29372. 
16295. 
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